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ABSTRACT 

We present 1.35 mm Submillimeter Array (SMA) observations around the 
low-mass Class source IRAS 19156+1906, at the the center of the LDN 723 
(L723) dark cloud. We detected emission from dust as well as emission from 
H2CO 3o,3-2o,2, DCN 3-2 and CN 2-1 lines, which arise from two cores, SMA 1 
and SMA 2, separated by 2'.'9 (880 AU in projected distance). SMA 2 is as- 
sociated with the previously detected source VLA 2. Weak SiO 5-4 emission 
is detected, possibly tracing a region of interaction between the dense envelope 
and the outflow. We modeled the dust and the H2CO emission from the two 
cores. The results from the modeling show that the cores have similar phys- 
ical properties (density and temperature distribution) but that SMA 2 has a 
larger P-H2CO abundance (by a factor 3-10) than SMA 1. The P-H2CO abun- 
dances found are compatible with the value of the outer part of the circumstellar 
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envelopes associated with Class sources. SMA 2 is harboring an active mul- 
tiple low-mass protostellar syst em and powering at least one molecular outflow 



(ICarrasco-Gonzalez et al.ll2008l ). In contrast, there are no known signs of out- 
flow activity towards SMA 1. This suggests that SMA 2 is more evolved than 
SMA 1. The kinematics of the two sources show marginal evidence of infall and 
rotation motions. The mass detected by the SMA observation, which trace scales 
of <1000 AU, is only a small fraction of the mass contained in the large scale 
molecular envelope, which suggests that L723 is still in a very early phase of star 
formation. Despite the apparent quiescent nature of the L723, fragmentation is 
occurring at the center of the cloud at different scales. Thus, at ~ 1000 AU the 
cloud has fragmented in two cores, SMA 1 and SMA 2. At the same time, at 
least one of these cores, SMA 2, has undergone additional fragmentation at scales 
of ^ 150 AU, forming a multiple stellar system. 

Subject headings: ISM: individual (LDN 723) — ISM: molecules — radio lines: 
ISM — stars: formation 



Introduction 



L723 is an isolated molecular cloud located at a distance of 3 00±150 pc (IGoldsmith et al. 



1984 ) and with a systemic velocity of Vlsr — 10.9 km s ^ ( Girart et al. 19971 ). It har- 



bors a low-m ass, Class 0, young ste llar object, first identified by the IRAS s atellite, IRAS 



19156-1-1906 ([Goldsmith et al.lll984j ). with a bolometric luminosity of 3.4 L© (jPartois et al. 



20051 ). The properties of the protostar in the IR to mm wavelength range have been ob 



ject o f several studies (lDavidsonlll987l : lAndre. et al.lll993l : iPezzuto et al.ll2002l : iDartois et al. 



20051). IRAS 19156-1-19 06 is associated with a CO outflow with a quadrupo l ar morphology 



1991 



(IGoldsmith et all Il984l : iMoriartv-Schieven k Snelll Il989l : lAverv et all Il990l : iHavashi et al. 



Lee et al.ll2002l ). This outflow consists of a pair of bipolar lobes aligned along the East- 
West direction (position angle PA ~ 100°) and another pair of bipolar lobes aligned roughly 
in the North- South direction {PA 32°), with IRAS 19156-1-1906 located at their common 



center ( lAvery et al.lll990l : iLee et al.ll2002l ). The East- West (EW) pair of lobes, and in par- 



ticular the blueshifted, eastern lobe, is a ssociated with several Herbig-Haro (HH) objects 
( iPalacios &: Eiroa]ll999l : iLopez et al.ll2007l ). Very Large Array (VL A) observations at 3.6 cm 
reveal two sources, VLA 1 and VLA 2, towar d the center of outflow (lAnglada et al.lll99lf) . al- 
though VLA 1 is likely a background source jAndada et al.lll996l : Ickart et al.lll997h . VLA 2 
shows the characteristics of a thermal radio-jet (jet-like morphology and partially optically 
thick free- free emission), and was first identified as the powering source of the EW pair of 
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molecular lobes (lAnglada et al.lll996l ). Recent, very sensitive, sub-arcsecond angular resolu- 
tion (0'.'2-0'.'7) VLA observations at 3.6 cm and 7 mm resolve VLA 2 into several components. 
The two brightest sources at 3.6 cm, VLA 2A and 2 B, are separated by 0^^3 (90 AU in projec- 



tion) and are possibly tracing embedded protostars (lAngladall2004j : ICarrasco-Gonzalez et al. 
20081 ). VLA 2 A is associated with extended emission along a PA ~ 115°. Two additional 



sources, believed to also trace embedded protos tars are identified at 7 mm: VLA 2C, lo cated 



Girart et al. 



at the position of the water maser emission ( 

VLA 2D, located ~ 3" southeast of VLA 2A. ICarrasco-Gonzalez et al. 



19971: iFuruva et al.l l2003D. and 
(|2008|) suggest that 



the CO high velocity emission and the HH objects are tracing three different outflows. One 
of them is the NS outflow possibly powered by VLA 2B. The other two are what previously 
was considered as the EW pair of lobes: one of the outflows has a PA ^ 115° and is powered 
by VLA 2 A, whereas the other one seems to be a "fossil" outflow with a PA ~ 90°. 



High angular resolution (3'.'5) VLA observations of L723 carried out by iGirart et al. 



(Il997l ). show that the NH3 emission arises from a V-shaped structure that traces the dense 
molecular envelope around the embedded protostars. This structure is elongated roughly 
east-west and is 0.15 pc long, with a mass of 7 Mq. The ammonia maps show evidence of 
heating and line broadening toward VLA 2. A second spot of heating i s observed 10 ^ ^ west 
of VLA 2. This western hot spot (WHS here after) is interpreted by iGirart et al.l (119971 ) 
as a very young protostar with and ag e shorter than 2 x 10^ yr. The NH3 structure is 
only partially traced by CS and N2H"'" (IHirano et al.l Il998l : IChen et al.l 120071 ). The dense 
envelope of L 72 3 wa s als o observed at sub-mil limeter wavelengths (850 fim and 450 /im) by 
Shirley et al.l ((2002) and lEstalella et al.l (120031 ) using the continuum SCUBA camera of the 
JCMT. A strong millimeter source was found, peaking at the position of VLA 2, with a large 
extension that matches the high-density gas traced by NH3. 

In this paper we present a study of the molecular and dust emission at scales of 
~ 1000 AU around the low-mass star forming region of L723. In § 2 we describe the 
observational procedure with the SMA. In § 3 we describe the main results obtained with 
the SMA data. In § 4 we present an analysis of the dust and H2CO emission. In § 5 we 
discuss the possible evolutionary scenarios of the sources in the region. 
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2. SMA Observations 

The observations were carried out with the 8-antenna SMA0 array at Mauna Kea 
at 1.3 millimeters in 2004 August (compact configuration). The phase center was set at 
a(J2000) = 19^17'^53!400 and 5(J2000) = 19°12'19'.'40. The phase calibrator used was QSO 
1925+211, and QSO 1749+096 was used as a reference calibrator. Absolute fiux and band- 
pass calibration was done by observing Uranus and Jupiter, respectively. The SMA correlator 
has a bandwidth of 2 GHz with 24 partially overlapping windows (chunks) of 104 MHz. All 
the chunks had 128 channels, except the second chunk which was configured to have 512 
channels. The correlator was setup to observe the H2CO 3o,3^2o,2 line in the second chunk 
of the lower side band. Maps were made with the (m, v) data weighted by the associated 
system temperatures and using a robust weighting of 0.0 (continuum data), 0.5 (H2CO) and 
1 (other line data). Table [1] lists the basic information about the resulting maps, including 
the frequency of the lines or continuum, the channel resolution (for line observations), the 
resulting synthesized beam and the rms noise of the maps. 



Results 



3.1. Dust continuum 



The 1.35 mm continuum map shows that the emission is clearly resolved into two com- 
ponents, L723 SMA 1 and SMA 2, with similar intensities (see Figure [1]) and separated by 
2'.'9 (880 AU in projected d istance) with a p osition angle of 106°. These two components 
were previously detected by iLaunhardtl (120041 ) at 3 mm. The strongest component, SMA 2, 
is associated with VLA 2 A, 2B and 2C, 3.6 cm and 7 mra continuum sources that are proba- 
bly young stellar object s (ICarrasco-Gonzalez et al.ll2008l ). The other component, SMA 1, is 
only detected at 3 mm (jLaunhardtll2004l ) and marginally at 7 mm (source VLA 2D: Carrasco 
et al. 2008) and does not have known counterparts at other wavelengths. Table [2] gives the 
values of position, intensity, fiux densities and the deconvolved sizes for both mm sources 
obtained by fitting a Gaussian to each component. The two sources are partially resolved at 
scales of about l'.'O-2'.'O (300-600 AU). 

Naturally weighted maps were obtained in order to better estimate the total flux of the 
double core and to try to detect the WHS. The total flux density detected by our observations 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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for SMA 2 and SMA 1 is 128 ± 6 mJy, which is 22 mJy higher (~20%) than that from the 
map shown in Figured], obtained with a robust weight of 0. This excess is possibly due to the 
better sensitivity of the natural weighted map to weak and extended emission. No emission 
is detected with the SMA at the position of WHS, with a 3-cr upper limit of 4.2 mJy beam~^. 
Gaussian tapers were applied to the visibilities to increase the beam size and see if the WHS 
could be detected, which yielded negative results. 

The mass of the circumstellar material traced by the dust emission at mm wavelength, 
where it is optically thin is: 

gS.D^ 



M 



1 



with g the gas-to-dust ratio, D the distance, 5*;^ the flux density, the dust mass opacity 
coefficient and By the Planck function for a blackbody of dust temperature Tdust. The 
major source of uncertainty in the mass derivation is t he poor know ledge of the distance 
(300 ± 150 pc). To derive the mass we adopted q = 10 (jDraind 120041 ). a dust mass opacity 
of ^250GHz = 0.9 cm^ g ^ (jOssenkopf fc Henning]|l994j ). To estimate t h e tern perature of the 
dust, we used the VLA NH3 (1,1) and (2,2) maps from iGirart et al.l (Il997l ). which have a 
slightly higher angular resolution (3'.'7 x 3'.'3) than the SMA maps. Following the standard 
analysis procedure for these two ammonia lines, we obtain that at the SMA 1, SMA 2 and 
WHS peak positions the ammonia rotation ten aperature is ^ 20 K , which converted to the 
gas kinetic temperature of the gas is ~ 25 K (jPanby et al.lll988l ). Table [3] gives the mass 
derived for the different components using this value for the dust temperature. The total 
mass of the emission detected by the SMA is ~ 0.2 4 M(^, which is o r ily a small fraction , 
~ 5%, of the dense envelope detected at larger scales (IGirart et al.lll997l : IShirley et al.ll2002l ). 
WHS has a very small mass, < 0.01 Mq. 



3.2. H2CO 

Figure [2] shows the channel maps for the H2CO 3o,3-2o,2 hne. The H2CO emission is 
concentrated around the two mm continuum sources, SMA 1 and SMA 2, although in most 
of the channels SMA 2 has brighter H2CO emission than SMA 1. Toward SMA 2 the H2CO 
emission peaks at or very close to the dust emission peak, whereas around SMA 1 the peak 
emission is, in some channels, displaced to the SE (only at wlsr = 10.79 km s~^ the H2CO 
peaks at the same position as the dust). The kinematics of the two sources is difficult to 
disentangle because their separation is similar to the angular resolution. The first order 
moment of the emission shows velocity gradients along the two cores with velocity shifts up 
to ~ 1 km s~^. Although there is not a clear simple trend, the emission east of SMA 1 
is clearly blueshifted with respect to the rest of the region. The velocity dispersion of the 
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H2CO emission (second-order moment map: see Fig. Hj) changes little around the two cores, 
with values within the range of 0.5 to 0.7 km s~^, which corresponds, for a Gaussian line, to 
a full width at half maximum of 1.2 to 1.7 km s~^. This velocity dispe rsion range is large r 



than the one measured from VLA NH3 observations, < 0.42 km s~^ (iGirart et al.l 119971 ). 
Figure [3] shows that the spectra of H2CO has larger line widths that those of NH3 (note that 
the (1,1) and (2,2) transitions have hyperfine components, which broadens the line, specially 
for the (1,1)). 



Since VLA 2 is a thermal radio jet with a well defined orientation (lAnglada et al.lll996l ). 
position-velocity (PV) plots along {PA = 116°) and perpendicular {PA = 26°) to the radio 
jet orientation were made crossing VLA 2, which lies at the peak of SMA 2 (Figure [5]). 
The PV plot along the radio jet direction passes very close to SMA 1 (see Fig. E]). In the 
direction perpendicular to the jet, the H2CO shows a velocity gradient (~ 0.8 km s~^ within 
1" around VLA 2) that could be indicative of rotation. In the direction along the radio jet 
there is overlap of the emission between SMA 2 and SMA 1. Nevertheless, at the position of 
SMA 1 the H2CO emission shows a relative minimum. In Figure [5] we also show the PV plot 
centered on SMA 1 with a position angle of PA = 30°, nearly perpendicular to the PV along 
the VLA 2 radio jet direction. This PV plot shows two weak peaks at different velocities, 
fLSPj^:^ 10.8 and 11.9 km s~^, that are slightly shifted in position (~ 0'.'9) . 

As in the case of the continuum emission, no emission is detected at the position of the 
ammonia WHS (see Fig. [3]). The 3-0" upper limit from the SMA map is 0.53 Jy beam~^ 
(2.4 K). 



3.3. SiO 

The SiO 5-4 line is marginally detected towards the SMA 1 and SMA 2 system (see 
Figure [6]). The emission appears clumpy and elongated in the NW-SE direction {PA ~ 145°), 
with the strongest emission about 11" SE of SMA 1. Taking into account the fact that the 
SiO is a tracer of molecular shocks, the emission appears relatively "quiescent": it is only 
slightly redshifted with respect to the systemic velocity (the emission is detected in the 
■^LSR = 11-31 and 12.44 km s~^ channels) and the line-width, Av ~ 2 km s~^, is only slightly 
larger than that of the lines tracing the dense circumstellar material such as H2CO and 
ammonia. 



The emission is appar ently more asso c iated with SMA 1, but given t he complex out- 
flow activity in the region (ILee et al.l 120021 : ICarrasco-Gonzalez et al.l |2008| ) it is difficult to 
elucidate with which outflow and powering source it is associated. Nevertheless, at the SiO 
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velocity and location there is strong CO emission (see channel maps from Lee et a. 2002), 
which suggests t hat the SiO is traci ng a region of interaction between the dense envelope 
and the outflow (ICodella et al.lll999l ). 



3.4. Other molecules 



Several hyperfine transitions of the CN 2-1 rotational transition were detected (Fig. [8]). 
The strongest line is a blend of three hyperfine lines, the F=7/2-5/2, 5/2-3/2 and 3/2-1/2 of 
the J=5/2-3/2 set. The J=3/2-5/2, F=l/2-3/2 hyperfine line is also marginally detected. 
Although there are other hyperfine transitions within the bandwidth their relative intensities 
were too small to be detected (see Fig. [8]). The integrated intensity of the detected transitions 
(Fig. [7]) shows that the emission is concentrated around VLA 2, being only weakly detected 
towards SMA 1. We also marginally detected the DON 3-2 line toward the two sources 
(Figs. Hand [HD- Despite the poor spectral resolution (~ 1 km s~^), these two molecules show 
emission at the ambient molecular gas velocity with narrow line widths (~ 1 km s~^). 



4. Analysis 

4.1. Model of the dust emission 

The dust emission from SMA 1 and SMA 2 are partially resolved with sizes of several 
hundred AU, which suggest that a significant contribution to the emission is possibly coming 
from the infalling envelopes around the disks. Figure [9] shows clearly that the correlated am- 
plitude in the visibility data centered on VLA 2 comes mostly from an extended component, 
although there is also a weaker unresolved component, with a flux density of ~ 20 mJy, that 
appears to be dominant at visibility radii larger than ~ 40 kA. To better characterize this ap- 
parently compact component, additional maps were obtained using only visibilities that have 
a radius in the {u, v) plane larger than 45 kA. The resulting map of the emission from SMA 1 
and SMA 2 appears to be compact. Both sources are unresolved with fluxes of 19.7 ±2.4 and 
13.3 ± 2.4 mJy for SMA 2 and SMA 1, respectively. This is about one third of the total flux 
density detected by the SMA for each source. A Gaussian flt to SMA 2 using the task IMFIT 
of AIPS yields an upper limit of the diameter of ~ 1'.'2: The other source, SMA 1, is too 
weak to do the same fltting. In any case, the compact dusty components seem to arise from a 
structure with a radius of < 180 AU. Without proper modeling it is not possible to elucidate 
whether this compact component traces the disk or the inner part of the envelope. Interest- 
ingly, the peaks of the two components are slightly shifted (~ 0'.'4) with respect to the peak 
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position given in Table [2l This could be produced by a somewhat asymmetrical distribution 
of the gas and dust in the envelope, which will not be a surprise because of the multiplicity 
detected with the VLA at the center of SMA 2. The co mpact component of SMA 2 coin cides 
well with the centimeter source VLA 2A detected by ICarrasco-Gonzalez et al. J2008h (the 
position of the dust is a(J2000)= 19'^17'^53^666 and 5(J2000)= 19°12'19'.'63). The peak of 
the SMA 1 compact component (a(J2000)= 19'^17'^53!920 and 5(J2000)= 19°12'18'.'52) is 
displaced to the northwest with respect to VLA 2D by 0'.'67 (see Fig. [TOl) . Given the weak 
detection of this source at both 1.35 and 7 mm (5 and A-a respectively) further more sensitive 
observations are needed to confirm this displacement. 

In order to characterize the properties of the double cores we modeled the dust emission 
assuming, for simplicity, that it arises from two independent optically-thin twin circumstellar 
envelopes with a radial density and temperature profiles at the position of the dusty sources 
detected with the SMA. The model integrates the dust emission assuming spherical symmetry 
from an inner radius, -Rmner, to an outer radius of 1.5 x 10^ AU (this value is somewhat 
arbitrary and does not affect the fit because the emission at these scales is filtered out by 
the SMA). We do not include the possible contribution from the circumstel lar disks, since it 



is po ssible that at this stage the disks are small, with a radius of < 30 AU, (IRodriguez et al. 



20051 ). The three possible pr otostars associated with SMA 2 are separated by a projected 
distance of ~ 0'.'9 or 270 AU (jCarrasco-Gonzalez et al.ll2008l ). below the angular resolution 
of our observations. If these sources are truly nearby protostars, it is possible that there is a 
cavity in the envelope between these sources. To check whether there is a significant cavity, 
we used different radii of 30, 100 and 180 AU. We adopted the density profile expected for 



an inf ailing envelope. 



n(r) oc r 



"1.5 



, up to an infall radius, -Rinfaii- Beyond -Rinfaii we used 



a r ^ density profile. In order to minimi ze the free pa r amete rs in the model, we adopted 



a -Rinfaii = 1000 AU, the value found by IShirley et al.l (120021 ). For a dusty cloud heated 
by an internal sourc e, the temperature radial profile can be characterized as T oc r~^/*^*^+^) 
(IKenyon et al.lll993l ). where /? is the dust emissivity spectral index. The value used is (3 = 1.5, 
which was derived from 850 and 450 fim SCUBA observations (lEstalella et al.ll2003l ). We 
adopted an 1.35 mm dust opacity of k ^ 0.9 cm^, which should be adequate for dusty 



particles with ice mantles at densities of about 10^ cm ^ according to lOssenkopf fc Henning 



(119941 ). Additional free parameters were the density nQ and temperature To at a radius of 
1,000 AU. Intensity profiles were obtained for each set of -Rmner, no and Tq and converted 
to a 2-D image map. These modeled maps were multiplied by the primary beam response 
of the SMA antennae, assumed to be 56" at 1.35 mm. Two sets of visibility data were 
generated: one obtained by subtracting the modeled map to the visibility continuum data 
of the SMA observations. The other set was obtained by replacing the values of observed 
visibilities with the values expected from the modeled map. Maps of the residual and of the 
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model were obtained from the visibilities using the same parameters as the SMA maps shown 
here. The rms and bias was finally computed for these residual maps. The bias is defined 
as the absolute fiux density of the residual on a given region (the region where the emission 
is detected). The models with an inner radius of 180 AU cannot reproduce well the SMA 
data, which suggests that if there is an inner cavity (produced by the presence of a multiple 
protostellar system), it is not big enough to affect the simple model used here. With respect 
to the solutions with an inner radius of 30 and 100 AU, there is, for both, a family of possible 
solutions in the Tg-rio plane. The plot of the rms for these two inner radii in the range of 
densities and temperatures computed are shown in Figure [TTl The best set of solutions can 
be expressed as uq = 1.5 x 10^ [To/30 K]-^-^ cm'^ and 1.1 x 10^ [To/30 K]-^ ^ cm^^ for the 
30 and 100 AU inner radius, respectively. 



The VLA ammonia maps from iGirart et al.l (119971 ) can be used to further constrain 



the model. At the positions of SMA 1 and SMA 2 the kinetic temperature estimated from 
the NH3 (1,1) and (2,2) maps is ~ 25 K (see § 13.11) . Since the beam of the VLA maps 
is ~ 3'.'5, it is reasonable to adopt the value of 1'.'75 or 520 AU, as the radius at which 
the gas is at 25 K. This implies that the temperature profile of the dust can be written as 
T = 20 (r/1000 AU)"*^'^^ K. The previous two equations that give the best set of solutions 
yield a density profile of ^(Hs) = 2.4 x 106(r/1000 AU)-^-^ cm'^ for i?in„er = 30 AU, and of 
1.8 X 10^(r/1000 AU)-^-^ cm'^ for i?inner = 100 AU. Figure [9] shows that the correlated fiux in 
the visibility domain for the model with -Rinner = 30 AU, Td = 20 K and no = 2.4 x 10^ cm~'^ 
fits reasonably well the SMA data. The synthetic maps from the model and from the residuals 
for this particular solution also show the remarkable similitude with the L723 dusty binary 
system (see Fig. [12]). Indeed, from the residual map obtained using only the longest baselines 
(see right panels of Fig. [T2!) the model of the binary envelope can account within the rms 
level of the map the compact, apparently unresolved, component of the two dusty sources, 
without the need of invoking the presence of an accretion disk (see the discussion section). 
On the other hand, the residual map done using all the visibilities shows some significant 
residuals (at a ~ 4-cr level): there are some positive residuals northeast of SMA 2 and 
extended negatives on the southern side of the dust emission from both SMA 1 and SMA 2. 
These residuals could be due to the departure from the spherical symmetry of the data. 



We can also compare the density we found with the value derived by iDartois et al. 



( I2OO5I ) by modeling the spectral energy distribution of L723. They estimated that at 100 
AU the density is 2.2 x 10"^ cm~^, which is a factor 2.6 to 3.4 lower than what we derive 
for To = 20 K. This difference could be due to the different dust opacity and the different 
density profile used. 

The derived density distribution can be used to estimate the mass of the two envelopes. 
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The equivalent radius of the SMA shortest basehnes in the visibihty domain (~ 12 kA) is 
3'.'75 (1125 AU). The combined mass for the two envelope within a radius of 1125 AU is in 
the 0.14-0.36 Mq range for a temperature range of Tq = 20-30 K and of -Rinncr = 30-100 AU. 
These values are compatible with the mass estimated in Table El 

As an additional test for the model, we have used a 3.2 mm continuum map obtained 
with the BIMA array. These observations were carried out in September 2003 in the C 
configuration. The continuum was observed simultaneously with the N2H"'" 1-0 and CH3OH 
2o-lo A and 2_i-l_i E fines (Masque, Girart & Estalella, in preparation). Tfie dust model 
was extrapolated to tfiis wavelengtfi, and tfie standard position of tfie BIMA antennas at 
tfie C configuration were used to create a syntfietic BIMA map using tfie same procedure as 
tfie one described for tfie SMA. Tfie syntfietic maps were obtained for tfie -Rjnner = 30 AU, 
Tq = 20 K and uq = 2.4 x 10^ cm~^ model. Figure [T^ sfiows tfiat tfie model also reproduces 
reasonably well tfie BIMA map. 



4.2. Modeling the H2CO emission 



Once tfie dust emission was success fully modeled, we modeled the Ho CO emission using 
the one dimensional version of Ratran ( IHogerheijde fc van der Takll2000l ). This is a Monte 
Carlo code that calculates the radiative transfer and excitation of molecular lines. Tfie code 
is formulated from tfie viewpoint of cells ratfier tfian pfiotons, wfiicfi allows tfie separation 
of local and external contributions of tfie radiation field. Tfiis give s an accurate and fast 
performance even for fiigfi opacitie s fjHogerfieijde &: van der Talal2000l ). Tfie H2CO collisional 
rates used were deriv ed by Green ( 199lh and were downloaded from tfie Leiden Atomic and 
Molecular Database ( Scfioier et al. 12005 ). 



We modeled tfie emission as arising from two cores witfi tfie same density and tem- 
perature profile, as suggested by tfie dust modeling. Since tfiere is a family of possible 
solutions, we fixed tfie temperature and density profiles derived in tfie previous section 
tfiat agrees witfi tfie ammonia rotational temperature: T = 20 (r/1000 AU)-^-^^ K and 



n(H2) = 2.4 X 10^(r/1000AU)-i-^ cm'^ (tfie solution for R 



in fall 



30 AU). 



Tfie kinematics is not well resolved witfi our observations, so we assume tfiat tfie gas is 
in free-fall collapse. Tfiis is marginally suggested by tfie line profile at tfie position of SMA 1 
and SMA 2 (see Fig. [T3] and § 15.11) . We adopted a velocity field for tfie free-falling gas of 
^^infaii = 0.5 (r/1000 AU)~*^'^ km s~^ and an intrinsic linewidtfi of 0.6 km s~^. Using tfiese 
values tfie modeled line widtfi was similar to tfie H2CO values measured witfi tfie SMA and, 
in any case, tfie use of tfiese values sfiould not affect critically tfie derivation of tfie H2CO 
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properties (|j0rgensen et al.ll2004l ) 



To run RATRAN we adopted a fixed H2CO abundance within a shell with an inner and 
outer radius and Rout, respectively, with the abundance being zero outside the shell. 
We ran a series of models with Ratran using a range of values for and Rout and for the 
P-H2CO abundance, X[p-H2C0]. Since the H2CO emission is significantly different between 
the two envelopes, we made synthesized maps using different sets of {Rm, Rout, X[p-H2C0]) 
values for SMA 1 and SMA 2. Synthesized maps with 16 channels and with a spectral 
resolution of 0.28 km s~^ were generated. We took into account that the systemic velocity 
of the H2CO between the two cores differs by ~ 0.2 km s~^ (see Table H]). The maps were 
multiplied by the primary beam of the SMA and then converted to visibilities. Maps of the 
model and of the residuals (generated by subtracting the model to the data in the visibility 
domain) were made in the same fashion as for the SMA dust emission maps (see § 14. ip . 



We explored values within the range of 8 x 10"^^ to 1 x 10"* for X[p-H2C0], 30 to 240 AU 
for Rin and 300 to 10^ AU for Rout- In order to find the best set of solutions, we computed 
separately the parameter for the spectra at the SMA 1 and SMA 2 positions and for the 
total flux density within a region of 11" x 8" centered around the twin envelopes. Table 
shows the range of the best fits found. The range of parameters is better constrained for 
SMA 2 than for SMA 1. This is due to its higher intensity and its apparently more compact 
size. Indeed, the model yields solutions that imply that the H2CO around SMA 2 should 
arise from the inner region of the envelope. More specifically, the outer radius cannot be 
significantly larger than 600 AU, but not lower than 300 AU (otherwise the emission would be 
unresolved, which is not the case). The SMA 2 abundance is also relatively well constrained, 
X[p-H2C0] ~ (3-10) X 10-1°. For SMA 1, the best set of solution yields H2CO abundances 
in the range X[p-H2C0] ~ (8-30) X 10"^^ but with values always lower than for SMA 2: 
the best solutions are found when SMA 2 abundances are a factor of 3 to 10 higher than for 
SMA 1. Another difference with respect to SMA 2 is that the best models for SMA 1 have 
a larger outer radius than for SMA 2, Rout — 600-5000 AU. 

Figure [15] shows the synthesized model, the SMA data and the residual for one of the 
best solutions found (see the figure caption for details). Figure shows the spectra toward 
SMA 1 and SMA 2 for the same model and for the SMA data. These figures show that the 
models fit the observations reasonably well, although the residuals appear up to about a 4-0" 
level. In spite of the good match between the synthesized and observed spectra for the two 
sources, the model reproduces better the morphology of SMA 2 than that of SMA 1. The 
morphology mismatch between SMA 1 and the synthesized map occurs for all the models. 
This suggest that the emission, specially from SMA 1, departs from the idealized spherical 
model used with the assumption that the H2CO arises only within a shell. Indeed, the 
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position-velocity cuts in two orthogonal directions centered on SMA 1 (see Fig. [5]) show a 
double peak in one of them and what roughly appear to partially be a "donut" , centered 
on SMA 1, in the other cut (partially because of the overlap with SMA 2 in the NW side). 
This suggests that the H2CO emission may not come from the shell of a spherical envelope 
but from a contracting and somewhat flattened shell or torus. For example, such a case was 
more clearly obser ved in the CS emission associated with the dense core ahead of HH SON 
(iGirart et al.ll200ll ). The marginal evidence that the H2CO arises from a torus supports the 
idea that the emission does not come from a spherical configuration. However, despite the 
limitations of the adopted approach, given the limited signal-to-noise ratio obtained and the 
use of only one transition, the results from the modeling should be regarded as a good first 
approximation to the properties of H2CO in these two cores. 



The derived P-H2CO abundances are compatible with the value found by |j0rgensen et al, 

(120051 ). 9.6 X 10"^*^, from single-dish observations, and are also within the range of values 



found fo r the so-called outer part of the c i rcurn stellar envelopes associated with Class 
sources (jj0rgensen et al.l l2005l : iMaret et al.l |2004J ) , which are defined as the region where 
temperatures are lower than ~ 100 K. In the inner region of the envelopes, where the 
temperature is > 100 K, the icy grain mantles evaporates and molecular abundances are 
expected to increase s i gnificantly: the H ^CO reaches abundances between 10~^ and 10~^ 
( lj0rgensen et al.ll2005l : iMaret et al.ll2004l ). Our model did not take into account this warm 
component. Nevertheless, the reasonable solutions found in this paper imply that the total 
amount of H2CO in the inner warm region is not significant, despite the very high H2CO 
abundance expected. This can be explained by the fact that for L723, the warm region is 
very small due to its low total bolometric luminosity, 3.4 L© (the temperature of 100 K is 
reached at a radius of 10 AU). 



5. Discussion 

5.1. The Twin Cores 

The dust emission reveals two similar substructures, with SMA 2 being only slightly 
more massive than SMA 1 (Table [2]). In fact, the modeling of the emission suggest that 
within the observational uncertainties (sensitivity and angular resolution) the dust properties 
of the two cores can be accounted as arising from "identical" twin cores. However, these 
twin cores have different molecular properties. The CN 2-1 and H2CO 3o3^2o2 emission is 
stronger in SMA 2 than in SMA 1. In fact, the modeling carried out in § 14.21 shows that 
the H2CO associated with SMA 2 has a higher molecular abundance (by a factor 3-10) and 
arises from a more compact region than in SMA 1. In both cores, the line width of the H2CO 
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lines are clearly larger (by up to 50%) than those of NH3. In addition, the H2CO 3o3-2o2 
emission is more compact than the NH3 emission, as can be clearly seen in Fig. [71 These two 
features taken together are a good indication of the presence of radial velocity gradients, with 
increasing velocity with decreasing radius to protostars. Both infall and rotation velocities 
are expected to have this behavior. The presence of a blueshifted peak slightky higher than 
the redshifted one at the position of the two cores (Fig. is an additional marginal evidence 
of infall motions. 

Previous observations find clear star formation signatures associated with the SMA 2 
core. First the presence of a well studied thermal radio jet, VLA 2 A, pinpoint the posi- 



tion of a protostar that is pow ering the east-west molecular outflow (lAnglada et al.lll996 



Carrasco-Gonzalez et al.ll2008l ). VLA high angular resolution observations reveal the pres- 
ence of a water maser, and 7 mm and 3.6 cm sources ( VLA 2B and 2C) ~ Gf!5 a round VLA 2A 



(IGirart et al.l 119971 : ICarrasco-Gonzalez et al.l |2008| ) . ICarrasco-Gonzalez et al.l (120081 ) inter- 
pret these as a young stellar system, which is forming in the 150 AU vicinity (in projected 
distance) around VLA 2A. In addition, they find that VLA 2B may be powering the north- 
south molecular outflow. This stellar system is being formed inside the SMA 2 core. In 
contrast, there are not known signs of outflow activity towards S MA 1, yet. There is only 



marg inal dust emission detected at 7 mm, labeled as VLA 2D (ICarrasco-Gonzalez et al. 
2008[ ). 



In conclusion, our study and analysis of SMA 1 and SMA 2 show that they have al- 
most identical properties (mass, and density and temperature distribution) and that there is 
marginal evidence of infall motions in both. However, the H2CO abundance towards SMA 2 
is signiflcantly higher than in SMA 1. Since SMA 2 shows clear outflow activity, in contrast 
with SMA 1, this suggests that SMA 2 is more evolved (due to a higher infall rate or because 
it is older). In any case, they are not "identical" twin cores but "fraternal" twin cores. 



5.2. The Mass Reservoir around the Class protostellar system 



The total mass derived from the SMA observations is less than 10% of the whole dense 



core, derived from SCUBA JCM T observations at 8 50 and 450 /im (lEstalella et al.ll2003l ) 
or from VLA NH3 observations (IGirart et al.l 119971 ). That is, the mass accreting to the 
protostars (or to the disks) at scales of < 1000 AU is only a small part of the dense core. 
In addition, from the modeling we found that the dust emission detected by the SMA at 
1.35 mm seems to arise from the dense twin envelopes, SMA 1 and SMA 2, and within the 
limits of our sensitivity there is no need to invoke to the presence of protostellar disks. That 
means that the contribution from the disks at 1.35 mm should not be higher than ~ 8 mJy 



-14- 



(the 5-cr value of the continuum residual map of Fig. [T2|) . Assuming that the disk has a 
temperature of 50 K, the upper limit of the disk mass is 6 x 10~'^ M©. This mass would be 
higher if the dust emission from the disk is not optically thin, but it would be lower if the 
temperature is higher. In any case, this suggests that the disk mass is less than 10% of the 
mass of each of the twin cores. Thus, the large scale (0.1 pc) molecular envelope traced by 
NH3 or SCUBA contain most of the mass or, in other words, it is still the major reservoir of 
mass to keep the star formation going on. Therefore, L723 seems to be in a very early phase 
of star formation. 

It is clear that fragmentation is occurring at different levels. The 0.1 pc scale molecular 
envelope has fragmented at the center and at scales of ~ 1000 AU into two cores, SMA 1 
and SMA 2. At the same time, at least one of these cores, SMA 2, has undergone additional 
fragmentation at scales of ~ 150 AU, forming a multiple stellar system. However, we note 
that within the uncertainties, the SMA 2 can be modeled a single envelope down to a 
radius of 30 AU. It is interesting to note that whatever is the fragmentation process that is 
acting in the different scales (turbulence, the initial magnetic field and angular momentum 
configuration) is occurring within a 0.1 pc dense envelope that appears to be v ery quiescent 



at sc ales of few thousands AU outside its center according to the NH3 data (iGirart et al. 



19971 ): line widths of only < 0.5 km s ^ and small velocity gradients, of 3 km s ^ pc 



6. Summary 

We present continuum and line 1.35 mm high angular resolution observations (~ 3") 
carried out with the SMA toward the low-mass protostellar region in the L723 dark cloud. 
We detected emission of the dust, the H2CO 3o,3-2o,2, DCN 3-2 and SiO 5-4 transitions, as 
well as a few lines of the CN 2-1 transition. We performed a radiative transfer analysis of the 
dust emission as well as for the H2CO line in order to constrain the physical and chemical 
conditions of the two cores. The main results of the paper can be summarized as: 

1. The dust emission arises from two similar cores, SMA 1 and SMA 2, separated by 2'.'9 
(880 AU in projected distance) and with masses of 0.06 and 0.08 M0, respectively. 
They are partially resolved, with deconvolved sizes of ~500 AU. SMA 2 is associated 
with VLA 2 and water maser emission. We modeled the dust emission from the two 
cores assuming that the cores have a density and temperature radial distribution of 
the type: n(r) cx r~^-^ and T cx r~^'^^. We found that within the uncertainty achieved, 
the two cores can be considered almost "identical" twins. The best set of solutions are 
those with a density and temperature that match the following expression: no(H2) = 
1.45 X 10^ [To/30 K]^'^ cm~^ {uq and Tq are the density and temperature at a radius 
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of 1000 AU). The model shows that all the emission detected can arise from the cores, 
with possibly no evidence of contribution from the accretion disk, which yields an 

0- 



upper limit of the disk mass of ~ 6 x 10 ^ M, 



2. The H2CO emission is concentrated around the two mm continuum sources, SMA 1 
and SMA 2, although the emission is brighter toward SMA 2. The kinematics of the 
two sources is difficult to disentangle because their separation is similar to the angular 
resolution. Nevertheless, there is marginal evidence of infall motions in both sources 
(double peak line, with brighter blueshifted component), of a velocity gradient around 
SMA 2 indicative of rotation, and that the H2CO in SMA 1 arises from a contracting 
and somewhat flattened structure. The linewidth of H2CO has values in the 1.2- 
1.7 km s~^ ran g e, wh ich is ~50% larger than the values measures for the NH3 lines 



by iGirart et al.l (119971 ). The H2CO emission was modeled adopting the density and 
temperature profile constrained by the dust emission. We used RATRAN assuming a 
fixed P-H2CO abundance within a shell with a inner and outer radius Rin and Rout, 
respectively, being the abundance zero outside the shell. Reasonable models were 
found for a P-H2CO abundance range of (3-10) x 10^^° and (8-30) x 10'^^ for SMA 2 
and SMA 1, respectively, but with the SMA 1 abundance lower than for SMA 2 by a 
factor 3 to 10. The best models are those that have a more compact emission in SMA 2 
{Rout ^ 300-600 AU) than in SMA 1 {Rout ^ 600-5000 AU). The P-H2CO abundances 
found are compatible with the value of the outer part of the circumstellar envelopes 
associated with Class sources. The total amount of H2CO in the inner warm region 
of the circumstellar envelopes (where T > 100 K) is not significant, despite the very 
high H2CO abundance expected (e.g.; Maret et al. 2004). This can be explained by 
the fact that for L723, the warm region is very small (a radius of < 10 AU) due to its 
low total bolometric luminosity, 3.4 Lq. 

3. The DCN and CN lines arise from the same region than the dust and H2CO. The 
CN emission peaks clearly toward SMA 2, whereas the DCN shows similar emission in 
SMA 1 and in SMA 2. 

4. SiO is detected marginally at slightly redshifted velocities around SMA 1 and SMA 2, 
and is possibly tracing a region of interaction between the dense envelope and the 
outflow. 



5. The Western Hot Spot (WHS), detected by lGirart et al.l (119971 ) from NH3 observations 



as a spot of local heating, is not detected by the SMA (line and continuum) observa- 
tions. The non-detection of the 1.35 mm emission yields an upper limit of the WHS 
mass of ~0.01 Mq. 
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6. The above results suggest that although SMA 1 and SMA 2 cores have almost identical 
physical properties (density and temperature distribution), they are rather "fraternal" 
twin cores, with SMA 2 being in a more evolved stage. Thus, SMA 2 is harboring an 
active multiple low-m ass protostellar system (VLA 2A , 2B and 2C), powering at least 
one molecular outflow (jCarrasco-Gonzalez et al.ll2008l ). In contrast, there are not signs 
of outflow activity towards SMA 1. 

7. The large scale, 0.1 pc, molecular envelope (traced by NH3 or SCUBA) is the major 
reservoir of mass to keep the star formation going on (it contains about > 90% of the 
mass) , which suggests that L723 is still in a very early phase of star formation. 

8. Fragmentation is occurring at different levels at the center of the L723 dense molecular 
cloud: at scales of ~ 1000 AU the cloud has fragmented into two cores, SMA 1 and 
SMA 2, and SMA 2 has undergon e additional fragmentation at s cales of ~ 150 AU, 
forming a multiple stellar system (jCarrasco-Gonzalez et al.l 120081 ). Whatever is the 
fragmentation process that is acting at the different scales (turbulence, the initial 
magnetic field and angular momentum configuration) it is occurring within a 0.1 pc 
dense envelope that appears to be at larger scales. 
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Fig. 1. — SMA map of the continuum emission at 1.35 mm (222.31 GHz) toward L723 
VLA 2. Contour levels are —3, 3, 5, 7, 9 , 11, 13, 15 and 17 times 1.8 mJy beam~^, the 
rms noise of the map. The synthesized beam is shown in the lower right corner. The crosses 
show the position of VLA 2A (the western cross) and VLA 2C, which are associated with 
SMA 2, and VLA 2D, associated with SMA 1 . These sources were detected by the VLA 
at 3.6 and 7 mm from sub-arcsecond maps by ICarrasco-Gonzalez et al.l (l2008l ). The titled 



cross shows the position of the WHS (iGirart et al. 
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Fig. 2. — SMA map of the H2CO 3o,3-2o,2 emission toward L723 VLA 2. Contour levels are 
—2, 3, 4, 5, 6, 7 and 8 times 0.18 Jy beam~^, the rms noise of the map. The synthesized 
beam is shown in the lower right corner of the lower right panel. The crosses show the 
position of the dust sources SMA 1 and SMA 2. The i'lsr velocity of the channels is shown 
in the lower left corner of the panels. 
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Fig. 3.— Spectra of the H2CO 3o,3-2o,2, NH3 (2,2) and (1,1) toward the position of SMA 2 
(left panels), SMA 1 (central panels) and WHS (right panels). The spectra were taken 
from channe l map s at similar angular resolution, ~ 3". The NH3 spectra are taken from 
Girart et al.l Jl997h . 
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Fig. 4. — Color image of the first order (bottom panel) and second order (top panel) 
moments of the H2CO emission towards L723 VLA 2 superposed with the zero order moment 
(integrated emission) contour map. The color image level is shown in the right side of 
the panels. The first contour is 0.1 Jy beam"^ km s~^, and the contour levels are 0.3 
Jy beam~^ km s~^. The crosses show the position of the dust sources SMA 1 and SMA 2. 
The synthesized beam of the maps is shown in the bottom right corner of the bottom panel. 



- 23 - 



13 



12 



11 



10 



— 1 1 1 r 

PA=26° 



-1 1 r- 




SMA 2 



—I (- 

13 I- PA=116° 



12 - 



H 1 1 1- 



11 - 



10 - 




— ' 

13 I- PA=30° 



^ 12 
E 



11 



10 



— 1 1 r— 




Offet position (") 



Fig. 5. — Position-velocity (PV) plots of the H2CO 3o,3-2o,2 emission. The top and central 
panels show the PV plot perpendicular and along the VLA 2 thermal radio jet direction, 
PA = 26° and 116°, respectively, centered on SMA 2 (and VLA 2). The bottom panel shows 
the PV plot centered on SMA 1 with a position angle, PA = 30° nearly perpendicular to the 
VLA 2 r adio jet and sim ilar to the north-south pair of lobes of the quadrupolar molecular 
outflow jLee et al.ll2002[ ). 
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Fig. 6. — Contour map of the SiO 5-4 emission overlaid with the grey scale image of the 
1.35 mm dust. Contours are —3, —2, 2, 3, 4 and 5 times 0.05 Jy beam~^, the rms noise of 
the map. The synthesized beam is sho wn in the bottom right corner. The filled triangles 
mark H2 knots ( jPalacios fc Eiroalll999l ). 
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Fig. 7. — Integrated emission of the H2CO 3o,3-2o,2, CN 2-1, DCN 3-2 lines observed with 
the SMA. For comparison we h ave also included th e integrated emission of the NH3 (2,2) 
line obtained with the VLA by iGirart et al.l (119971 ). The synthesized beams are shown in 
the lower right corner of each panel. Contours are —2, 2, 3, 4, 5, 6, 7, 8, 10 and 12 times 
3.2, 34, 16 and 51 mJy beam'^ km s'^ for NH3, DCN, CN and H2CO, respectively. The 
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Fig. 8.— Spectra of the H2CO 3o,3-2o,2, DCN 3-2, CN N =2-1, J = 5/2-3/2 and J = 3/2- 
1/2 obtained averaging a region of 7" x 5" around SMA 1 and SMA 2. The CN and DCN 
spectra has a lower velocity resolution (~ 1.0 km s~^). For the CN hyperfine transitions, 
a vertical dashed line indicates the position of the expected hyperfine transition with the 
height proportional to its relative intensity. 
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Fig. 9. — SMA correlated flux for the 1.35 mm continuum emission versus uv distance (in 
units of kA) at the position of SMA 2. The correlated flux was derived by vector averaging 
the amplitude of the visibilities over annular bins. The bins have a width of 3.3 kA for radius 
lower than 50 kA and of 10 kA for larger radii. The solid line shows the expected value for 
the amplitude assuming no signal (i.e., the "zero bias"). The dotted line shows the expected 
flux for an unresolved source of 20 mJy. The open squares show the expected correlated for 
the model Tq = 20 K and no(H2) = 2.4 x 10*^ cm"^ described in § SJl 
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Fig. 10. — SMA map of the 1.35 mm dust emission obtained using only the visibihties 
with a radius larger than 45 kA in the {u, v) plane and applying a Natural weighting to the 
visibihties. The contour levels are —3, 3, 4, 5, 6 and 7 times the rms noise of the map, 
2.4 mJy beam~^ The synthesized beam, 2'.'94 x 0'.'92 with a PA = 74.2°, is shown in the 
bottom right corner. The crosses show the position of VLA 2A (the western cross) and 
VLA 2C, which are associated with SMA 2, and VLA 2D, associated with SMA 1. 
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Fig. 11. — Plots shows the locus of best solutions represented by the function of the rms of 
the residual map derived by subtracting the dust model for an inner radius of 30 {left panel) 
an 100 AU {right panel) to the SMA visibility data of the 1.3 mm continuum emission. Levels 
are in steps of 3 times the rms of the original SMA map, 1.5 mJy beam"^. The dashed 
grey line shows the line of approximately set of best solutions represented as the function 
n(H2) = 1.5 X 106[T/30K]^-2 cm'^ (left panel) and ^(Ha) = 1.1 x 10^ [r/30 K]^-2 cm'^ 
(right panel). 
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Fig. 12. — SMA synthetics maps for the dust model 



30 AU, Tn 



20 K and 



no = 2.4 X 10^ cm~^ described in § 14.11 Left panels: the model (top) and residual (bottom) 
dust maps obtained using all the visibilities observed with the SMA and using robust of 1. 
Contours levels for the model map are are —3, 3, 5, 7, • ■ • ,21 times 1.5 mJy beam~^, 
the rms of the map. Contour levels for the residual map are —3, —2, 2, 3, and 4 times 
1.5 mJy beam~^, Right panels: the model (top) and residual (bottom) dust maps obtained 
using visibilities with a radius larger than 45 kA and Natural weighting (as in Fig. [TOj) . 
Contour levels are —3, —2, 2, 3, 4, 5, 6 and 7 times 2.4 mJy beam~^. The synthesized 
beams are shown in the lower left corner of each panel. The crosses show the position of 
SMA 1 and SMA 2. 
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Fig. 13. — Left panel: BIMA map of the 3.2 mm continuum emission obtained with the C 
configuration in September 2003. The synthesized beam, 10'.'6 x 6'.'7 and PA = —7°, is shown 
in the bottom right corner. Contours are —2, 2, 3, 4, 5, 6 and 7 times 2.8 mJy beam~^, 
the rms noise of the map. The crosses show the position of SMA 1 and SMA 2. Middle 
panel: BIMA synthetic map at 3.2 mm continuum at the C configuration for the dust model 
To = 20 K and n(H2)o = 2.4 x 10^ cm~^ described in § 14.11 Contours are the same as in the 
previous panel. Right panels: BIMA residual map obtained from the previous two maps. 
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Fig. 14. — Synthetic (solid line) and observed (dashed grey line) H2CO 3o,3-2o,2 spectra 
at SMA 1 (left panel) and SMA 2 (right panel). The model parameters for the synthetic 
spectra are the same as in Figure [T5l 
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Fig. 15. — Channel maps of one of the best models, the SMA data and the residual 
(SMA— model) for the H2CO 3o,3-2o,2 emission. Contours are —3, —2, 2, 3, 4, 5, 6, 7 
and 8 x 0.17 mJy beam"^. The synthesized beam is shown in the lower right corner of the 
lower right panel. The crosses show the position of the dust sources SMA 1 and SMA 2. The 
parameters for this model are ^[p-HaCO] = 3 x 10~^, = 240 AU and Rout = 600 AU for 
SMA 2, and ^[p-HaCO] = 1 x 10-^°, = 60 AU and Rout = 5000 AU for SMA 1. The 
^LSR velocity of the channels is shown in the lower left corner of the right panels. 
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Table 1. Parameters of the SMA observations 







Synthesized Beam 




Spectral 


rms 




V 


HPBW 


PA 


Az/ 


Resolution 


Noise 


Observation 


(GHz) 


(arcsec) 


(deg) 


(MHz) 


(km s~^) 


(mJy beam"-*^) 


Continuum 


222.31 


3.3 X 1.5 


79 


4000 




1.8 


H2CO 3o,3-2o,2 


218.2222 


3.1 X 1.8 


76 


104 


0.28 


175 


DCN 3-2 


217.2385 


3.2 X 2.1 


76 


104 


1.07 


78 


SiO 5-4 


217.1050 


3.2 X 2.1 


76 


104 


1.12 


60 


CN 2-1 5/2-3/2 7/2-5/2^^ 


226.8748 


3.1 X 2.0 


76 


104 


1.07 


65 


CN 2-1 3/2-1/2 5/2-3/2^^ 


226.6595 


3.1 X 2.0 


76 


104 


1.07 


88 



^CN N,J,F hyperfine transitions 



Table 2. Dust properties at 1.35 mm^ 





Positions 






Deconvolved Size 


PA 


Source 


q;(J2000) 5(J2000) 


(mJy beam~^) 


(mJy) 


(arcsec) 


(deg) 


SMA 1 
SMA 2 


19''17"53!884 19°12'18'.'51 
19''17'"53!694 19°12'19'.'68 


29.9±1.8 
33.7±1.8 


45±4 
59±4 


(2.0 ±0.3) X (<0.6) 
(2.3 ±0.2) X (1.0 ±0.2) 


56 ±9 
73 ±6 



^Values corrected by the primary beam of the SMA antennas. 
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Table 3: Masses of the different components in L723 





5^(f.35mm) 


^dust 


M 


Source 


(mJy) 


(K) 


(Me) 


All" 


128±6 


25 


0.24 


SMA 1 


45±4 


25 


0.09 


SMA 2 


59±4 


25 


0.11 


WHS 


< 5.7 


25 


< 0.011 



" To estimate the total mass we took into account the contribution of the mass from the 
24 mJy excess of emission of the Natural map (obtained adopting a temperature of 
20 K). 



Table 4: H2CO 3o,3-2o,2 









^LSR 


Av 


Source 


(K) 


(K km s^i) 


(km s^^) 


(km s~^) 


SMA 1 


4.0±0.7 


8.7±1.0 


11.47±0.11 


2.0±0.3 


SMA 2 


8.7±0.7 


14.6±0.8 


11.64±0.04 


1.6±0.1 



Table 5: Best H2CO Models 









Rout 


Source 


X[p-H2C0] 


(AU) 


(AU) 


SMA 1 


(8-30) X 10~" 


30-120 


600-5000 


SMA 2 


(3-10) X 10-1° 


30-240 


300-600 



